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ABSTRAGT 

I propose that some of the most luminous planetary nebulae (PNs) are actually proto-PNs, where a 
companion white dwarf (WD) accretes mass at a relatively high rate from the post-asymptotic giant 
branch star that blew the nebula. The WD sustains a continuous nuclear burning and ionizes the 
nebula. The WD is luminous enough to make the dense nebula luminous in the [O HI] A5007 A line, In 
young stellar populations these WD accreting systems account for a small fraction of [O HI]-luminous 
PNs, but in old stellar populations these binaries might account for most, or even all, of the [O HI]- 
luminous PNs. This might explain the puzzling constant cutoff (maximum) [O HI] A5007 luminosity 
of the planetary nebula luminosity function across different galaxy types. 

Subject headings: binaries: close—planetary nebulae: general—stars : AGB and post-AGB 


1. INTRODUCTION 

The planetary nebula (PN) luminosity function 
(PNLF) depends on the age of the parent population 
(Dopita et al. 1992; Mendez et al. 1993; Marigo et al. 
2004), and to lesser degree on the metallicity (Giardullo 
& Jacoby 1992; Marigo et al. 2004), hence it varies be¬ 
tween different types of galaxies (e.g., Giardullo et al. 
2004). However, the [O HI] A5007 A (hereafter [O HI]) 
most luminous end of the PNLF, and in particular the 
cutoff (maximum) [O HI] luminosity, seems to be similar 
in all large PN populations, with very small dependance 
on galaxy type (Giardullo et al. 2005). This allows a very 
successful use of the PNLF as a standard candle, from 
the galactic bulge (Pottasch 1990), through the LMG and 
SMC (Jacoby et al. 1990), M31 (Giardullo et al. 1989) 
and to galaxies at larger distances, spirals (e.g. Feldmeier 
et al. 1997), and ellipticals (Jacoby et al. 1996). 

Some basic properties of the PNLF are well understood 
(e.g., Dopita et a. 1992; Mendez et al. 1993; Marigo et 
al. 2004; review by Giardullo 2006). However, a major 
puzzle is the cutoff (maximum) luminosity in old stellar 
populations, such as in elliptical galaxies, which is the 
same as in young populations (Giardullo et al. 2005; 
Giardullo 2006). This cutoff luminosity of the PNLF is 
L[0 HI] ~ 6OOL0, and it leads to a chain of constraints 
(Giardullo et al. 2005; Giardullo 2006): Such [O HI] 
luminosity requires the ionizing central star to have a 
bolometric luminosity of L* > 6000 Lq, which in turn 
requires the central star to have a mass of M* > O.QMq, 
and the progenitor to have a main sequence mass of M > 
2Mq. Such a progenitor mass is not expected in old 
stellar populations. Single star evolution alone cannot 
account for this finding (Marigo et al. 2004). 

Giardullo et al. (2005) propose that the most [O HI] 
luminous PNs in old populations are descendant of blue- 
straggler type stars; namely, two lower mass stars, ~ 
IM0, merged on the main sequence to form a star of 
mass ~ 2Mq. Such a population of luminous PNs is 
quite possible. I propose another plausible type of [O HI]- 
luminous PNs, composed of a post-AGB star and a close 
accreting white dwarf (WD) companion. 


2. OBSERVATIONAL HINTS 

Several types of observations hint at the possibility 
that some of the [O HI]-luminous PNs are actually proto- 
PN systems ionized by an accreting WD companion. 

(1) Morphology: Throughout this paper it is assumed 
that the bipolar morphology of PNs is acquired via bi¬ 
nary interaction (Soker 1998). I stress that the stellar 
binary model for shaping bipolar PNs is a conjecture, 
and that it is not accepted by all researchers in the field. 
In the local galactic PNLF according to Mendez et al. 
(1993) four of the five most luminous PNs are bipolar, 
NGG 7027 (Latter et al. 2000), NGG 7026 (Gorradi & 
Schwarz 1995); NGG 6572 (Miranda et al. 1999) and 
NGG 6537 (Gorradi & Schwarz 1995). Only 7009, fourth 
most luminous, is elliptical but with two jets along the 
symmetry axis (Balick 1987). We note that NGG 7027, 
the most luminous by a large margin in the Mendez’s 
sample, is a near by compact PN, only recently turned 
into a bipolar PN (Kastner 2005). At larger distances, 
several kpc or more, its bipolar structure would have 
been unresolved, and it would have been classified as a 
round or an elliptical PN. Therefore, the fact that many 
[O HI]-luminous PNs in the Magellanic clouds are classi¬ 
fied round or elliptical (Magrini et al 2004; Stanghellini 
2005) is not surprising, and cannot be used as a contra 
observation to the claim above. Some spherical [O HI]- 
luminous PNs can be formed from massive single stars. 
However, I expect many of the [O HI]-luminous PNs that 
are classified round and elliptical to have inner bipolar 
structure, as in NGG 7027. Blue stragglers are not ex¬ 
pected to form bipolar PNs (Soker 1998) unless a close 
tertiary star exists, but this is very rare. It should be 
noted that in the galaxy and the Magellanic clouds there 
are young stellar populations, and most of the [O HI]- 
luminous PNs can be explained by single star evolution 
(Marigo et al. 2004). The point here is that I expect a 
large fraction of the most [O HI]-luminous PNs in these 
galaxies to have a binary companion, with some of the 
companions being accreting WDs. 

(2) Young PNs: The central star of NGG 7027 left the 
AGB ^ 700 years ago (Latter et al. 2000). The bipolar 
PN NGG 7026 is at a similar evolutionary stage (Cuesta 
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et al. 1996) with a dynamic age of ~ 800 years (Solf 
& Weinberger 1984). The age of the luminous bipolar 
PN NGC 6572 is similar to that of these two luminous 
PNs (Miranda et al. 1999). These young PNs suggest 
that many [O III]-luminous PNs are young. Interestingly, 
NGC 7027 has a relatively high extinction (Robberto et 
al. 1993; Kastner et al. 2002), showing that a PN does 
not need to clear all its surrounding dust before becom¬ 
ing a luminous PN. In the Milky Way galaxy, unlike in 
elliptical galaxies with old stellar populations, there are 
massive progenitors that can explain the most [O III]- 
luminous PNs. However, even in these a fast evolution is 
required for the PN to spend a long enough time period in 
an ionized bounded phase. In the ionized bounded phase 
all ionizing radiation is absorbed by the nebula; hence, 
the formation of the [O III] line is more efficient. A bi¬ 
nary companion can disturb the AGB envelope, causing 
high mass loss rate and fast post-AGB transition to the 
PN stage. 

(3) Luminous supersoft X-ray sources (SSS): Many of 
the SSS are WDs accreting mass and sustaining continu¬ 
ous nuclear burning. A large fraction of these WDs have 
a luminosity of Two ~ lO^L©, and a temperature of 
Twd ~ I — 2 X 10® K (e.g., Di Dtefano & Kong 2004), 
appropriate for the formation of [O III]-luminous PNs 
(Dopita et al. 1992). We note that in the WD-accreting 
model proposed in this paper the WDs accrete at higher 
rates than most SSSs. However, because the soft X-rays 
of SSSs are readily absorbed by material shed by the 
post-AGB primaries, SSSs in close binary systems with 
post-AGB primaries are likely to go undetected in X- 
rays. It should be noted that in the present model close 
means orbital separation of ~ I AU, and not ~ few x Rq. 
Hence, the binary systems that will formed from such 
systems will not merge to form type la SNs. 

(4) Central close binary systems in PNs: Some binary 
post-AGB stars demonstrate that binary systems at an 
orbital separation of ~ I AU, with main sequence or WD 
companions, can survive without entering a common en¬ 
velope phase (Van Winckel 2004). The companion leads 
to the formation of a bipolar nebula, as in the Red Rect¬ 
angle (Gohen et al. 2004). The presence of a disk around 
the Red Rectangle and similar systems suggests that this 
phase can last hundreds of years and longer. 

(5) Symbiotic systems evolving to PNs: There are sev¬ 
eral similarities between symbiotic nebulae and bipolar 
PNs (Gorradi & Schwarz 1995; Kwok 2003), causing 
some misclassifications. However, many of the symbiotic 
systems are evolving to become PNs. In a PN, by defini¬ 
tion, the star that ionizes the nebula is the star that blew 
the nebula. In symbiotic nebulae the secondary ionizes 
the nebula as they accrete mass from the primary mass¬ 
losing giant star. Schwarz & Gorradi (1992) suggested 
that the bipolar nebula BI Crucis is an object evolv¬ 
ing to become a PN, where there is a central accreting 
WD. The type of objects discussed here are similar to 
their model of BI Grucis. In M2-9 the primary star 
is most probably a post-AGB star (Lim & Kwok 2000), 
while the nebula is ionized by an accreting secondary WD 
(Goodrich 1991; Schwarz et al. 1997). Glose examination 
suggests that M2-9 is a symbiotic nova type system, and 
not a PN (Solf 2000). However, had it been located at 
a farther distance, such as the galactic bulge, it would 
have most probably been classified as a PN. This classi¬ 


fication would be even more probable was M2-9 located 
in external galaxies. As the mass losing star begins to 
ionize the nebula, this system most likely will go through 
a smooth transition to a bona-fide PN. 

Based on the observations listed above, I propose that 
some of the [O HI]-luminous PNs are binary systems 
where one of the stars, the primary, is a post-AGB star, 
while the secondary star is a WD accreting from the pri¬ 
mary and sustaining a continuous nuclear burning. Such 
a system is actually a symbiotic system, because the ion¬ 
izing source is the companion to the mass-losing star. 
However, at this stage of evolution and at the large dis¬ 
tances of the galactic bulge and external galaxies, the 
dense ionized nebula has all the characteristics of PNs 
(Kwok 2003). As suggested by the observations of near 
by proto-PN systems, this phase can last for hundreds of 
years. These systems are similar to M2-9 but the orbital 
separation is ~ I AU, like that in the Red Rectangle, 
rather than ~ 27 AU as in M2-9 (Livio & Soker 2001). 
The much smaller orbital separation than that in M2-9 is 
required for a much higher accretion rate by the WD (sec. 
3.2). The high mass accretion rate is required to sustain 
a continuous nuclear burning, causing the system to be¬ 
come an SSS. However, because of the relatively high 
accretion rate, the WD swells and becomes somewhat 
cooler, and because the post-AGB star loses mass at a 
higher rate than the mass losing star in most SSSs, the 
SSS is highly obscured. In external galaxies these types 
of systems most likely will not be classified as SSSs. 

3. MAIN PROPERTIES OF THE PROPOSED MODEL 
3.1. Relatively high mass accretion rates 

To sustain a continuous nuclear burning the WD ac¬ 
cretion rate should be higher than some critical value 
which increases with WD mass; this critical mass ac¬ 
cretion rate is Mwdc — 8 x yr“®, for a WD 

mass of Mwd = 0-QMq, Mwdc — 1-6 x 10~'^Mq yr“^ 
for Mwd = O.SMq, and Mwbc — 2.5 x yr“®, 

for Mwd = lAf© (Hachisu et al. 1999). At an accre¬ 
tion rate of Mwd ~ (2 — 4)Mwdc the accreting WD 
swells, blows a wind, and its temperature drops from 
PwD > 3 X 10® K to lower temperatures (Hachisu et al. 
1999). For example, a temperature of Twd = 2 x 10® K 
(this is the temperature of the central star of the PN 
NGG 7027; Latter et al. 2000) is reached at accretion 
rates of Mwd — 1-6 x IO^^Mq yr“^ for Mwd = 0.6Mq, 
Mwd — 6 X IO^’^Mq yr“^ for Mwd = O. 8 M 0 , and 
Mwd — 1 O“®M 0 yr“^ for for Mwd = IM 0 . As tem¬ 
perature drops toward Twd ~ 1.2 x 10® K the formation 
of the [O HI] line becomes more efficient (Dopita et al. 
1992). 

3.2. Close binary systems 

For the WD to accrete at a high enough rate, the or¬ 
bital separation cannot be larger than a few AU. Gon- 
sider the post-AGB phase and scale the central star post- 
AGB stellar wind parameters following the properties 
of M2-9. The primary mass loss rate drops already to 
Magb ~ few X 1 O“®M 0 yr“^, but the post-AGB enve¬ 
lope is still large with a wind speed of ~ 20 km s“^. 
With the orbital speed of ~ 30 km s“^, the WD-wind 
relative speed is Urei ~ 35 km s“^. The WD accretion 
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rate from the wind is 
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This implies that the WD accretes a large fraction of the 
post-AGB wind. The luminosity is Two > lO'^T©, with 
part of the accreted mass being blown by the WD, and 
the WD swells to a larger radius (Hachisu et al. 1999). 
For a given accreting WD sustaining a continuous nuclear 
burning the luminosity is constant at Twd = 5000L© for 
lUwD = O.GMq, = lO'^T© for = 0.7Mq, and 

UwD = 3 X lO'^L© for Mwd = IIU© (van den Heuvel et 
al. 1992). 

Considering the critical mass accretion rate given in 
§3.1, equation m shows that the orbital separation must 
be a < 1 AU. Such binary systems must experience 
strong binary interaction, namely, strong tidal interac¬ 
tion, and/or Roche lobe overflow (RLOF), and/or a com¬ 
mon envelope interaction. 


3.3. Bipolar morphology 

For the [O III] line not to be obscured by dust, a col¬ 
limated fast wind (CFW) or jets blown by the accreting 
WD should clear a path with lower optical depth, as in 
the model for M2-9 proposed by Livio & Sober (2001). 
As evident from the case of NGC 7027 (§2) it is enough 
that extinction is low in only part of the nebula. Indeed, 
as discussed in connection to M2-9, the accreted mass 
has a high enough specific angular momentum and an 
accretion disk is formed around the WD. Such a disk is 
assumed to blow jets. In the case of a common envelope 
interaction, the formation of a bipolar nebula requires 
the WD to blow jets before or after the common envelope 
phase (Sober 2004). An accretion disk around the WD 
will be formed via RLOF. An accretion disk is likely to 
be formed, and blow jets, even when accreting from the 
AGB or post-AGB wind. The condition for the formation 
of an accretion disk around the WD can be expressed in 
the form Rd > i?wD, where Rd is the outer radius of the 
disk that would form given the wind’s specific angular 
momentum. The disk radius is given approximately by 
(e.g. Livio & Warner 1984; see Livio & Sober 2001 for 
the full expression) 


Rd ~ Rq 



V M© j 


’^rel 


50 km s 


-1 


cm. 


( 2 ) 

An Mwd = 0.6M© WD swells to a radius of i?wD ^ 177© 
as long as the accretion rate is Mwd 2 x 10“®M© yr“^. 
For higher mass WDs the mass accretion rate should be 
even higher. We conclude that for the conditions as¬ 
sumed here, where the mass loss rate from the primary 
already dropped from its maximum value on the upper 
AGB, an accretion disk will be formed, leading to jets 
(or a GFW) and the formation of a bipolar PN. 


3.4. Mass loss rates 

Still, for the GFW from the WD to clear a path 
for the ionizing radiation and ionizing a large volume 
of the nebula, the primary post-AGB mass loss rate 


cannot be too high; as seen below, the constraint is 
Magb ^ 10“®M© yr“^. At the same time the mass loss 
rate should be Magb ^ 10“®M© yr“^, depending on the 
orbital separation and WD mass, to sustain continuous 
nuclear burning on the accreting WD. 

Following the analysis of Livio & Sober (2001) of M2-9, 
I compare the nebular recombination rate with ionizing 
photon emission rate. The total recombination rate from 
a minimum distance rmin to large distances in the nebula 
ri Cinin is given by (assuming a spherical flow) 


Magb 


R = A-k 
, 2 


10“® M© yr’ 




20 km s 


ariiTier'^dr ~ 5 x 10^^ 


where a is the recombination coefficient, rie and rii are 
the electron and ion number densities (respectively). 
Here we assumed a fully ionized nebula of solar compo¬ 
sition. The number of ionizing photons emitted by the 
WD per unit time, assuming a temperature of Twd — 
2 X 10’^ K, is 




( 4 ) 


The two rates obtained above, for R and S, are within 
an order of magnitude from each other, implying that 
the WD can ionize the nebula to moderate distances in 
cases of a spherical geometry. Jets (or a GFW) which are 
expected to be blown in the WD-accreting model, can 
clear a paths through the dense slow wind to a distance 
of ~ 3 — 10 times the orbital separation. Namely, rmin ~ 
3 — 10 AU, hence the recombination rate is much lower 
than that given by equation In this bipolar geometry 
the nebula can be ionized to larger distances than in the 
spherical geometry, even for mass loss rate as high as 
Magb — I0“^M© yr“^. The extended nebula regions 
are less obscure by dust, enabling the escape of the [O III] 
line. 


3.5. Rapid post-AGB evolution 

As mentioned in §3.2 the constraint of high accretion 
rate from a star after its AGB peak mass loss rate implies 
a close binary system. Such strongly interacting WD 
companions cause large disturbances in the post-AGB 
envelope of the primary star, leading to an enhanced 
mass loss rate and a fast transition from the AGB to the 
PN phase. However, while the WD is accreting at a high 
rate, it is the WD secondary that ionizes the nebula, not 
the primary. Therefore, by definition this is a symbiotic 
like system and not a PN. From the observational point 
of view, it might be classified as a young PN, as M2-9 
is. Because of the fast evolution and small orbital sepa¬ 
ration, the primary will start ionizing the nebula shortly 
after the WD ceased to ionize it, or even before the WD 
ceased to ionize the nebula. When the primary star ends 
with a low mass core, such as when the progenitor is of 
low initial mass, the most luminous phase occur when 
the accreting WD ionizes the nebula, and not the PN 
phase itself. 


3.6. Statistics. 
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Ciardullo et al. (2005) suggested that the progenitor 
of most [O III]-luminous PNs in old stellar populations 
are blue stragglers, with a mass about twice the main- 
sequence turnoff mass. They estimate that these PNs 
spend Tb — 500 yr in the luminous stage. With this value 
of Tb they estimate that ~ 5% of stars should become 
[O III]-luminous PNs. They then suggest that the frac¬ 
tion of blue stragglers in old-population galaxies is well 
suited to the proposition that all [O III]-luminous PNs 
in these galaxies result from blue stragglers. In young 
stellar population galaxies, such as the Milky Way, most 
[O III]-luminous PNs comes from massive stellar progen¬ 
itors. One of the problems I see in Ciardullo et al. (2005) 
estimation is that most blue stragglers have masses be¬ 
low twice the main sequence turnoff mass (e.g., Lee et 
al. 2003; Piotto et al. 2004; De Marco et al. 2005), 
with less than ~ 10% having a mass more than 1.7 times 
the turnoff mass. The final core mass of the PN de¬ 
scendants of these blue stragglers, though, will be higher 
than that of main sequence star of the same mass. It is 
not clear that there are enough blue stragglers to explain 
the [O III]-luminous PNs in old-population galaxies. 

Soker & Rappaport (2000) estimated that in the Milky 
Way galaxy ~ 1% of PNs have been formed in a binary 
system where the secondary was a strongly interacting 
WD, as required for the WD-accreting model proposed 
in the present paper. This percentage could be some¬ 
what larger, because Soker & Rappaport did not evolve 
systems which had strong interaction on the first giant 
branch (RGB); some of these might eventually also form 
PNs progenitors with strongly interacting WD secondary 
star. Also, in old stellar populations the fraction of WD 
companions is expected to increase, as the massive com¬ 
ponents of binary systems have already gone through the 
PN phase. Therefore, it seems that PNs progenitors 
with strongly interacting WD companions can account 
for ~ 1 — 3% of the PN population. The fraction of the 
systems discussed here among all observed PNs can be 
even higher if the vast majority (Soker & Subag 2005), 
or even all (De Marco & Moe 2006), PNs are formed by 
binary stars. 

In addition, the accreting WD companion can stay lu¬ 
minous for longer than the Tb = 500 yr assumed by Cia¬ 
rdullo et al. (2005). As an example consider again M2- 
9. The accreting WD companion is currently still ion¬ 
izing the nebula (Goodrich 1991; Schwarz et al. 1997). 
The dynamical age of the outer nebula is > 1200 yr 
(Schwarz et al. 1997; Solf 2000), suggesting the the bi¬ 
nary interaction, including high mass transfer rate, can 
last for > 1000 yrs. Schwarz et al. (1989) suggest that 
He 2-104 is both a proto-PN and a symbiotic system. 
They estimate the nebular dynamical age to be 800 yrs. 
These studies of M2-9 and He 21-104 show that a sys¬ 
tem can be a symbiotic type, with a high mass loss rate 
from the primary simultaneously with ionization by an 
accreting WD companion, for > 1000 yrs. This time 
scale is compatible with that expected in the model. A 
post-AGB star shrinks to ~ 50% its maximum AGB size 
when its envelope mass is ~ 3 x IO^^Mq, while the rapid 
shrinking occurs when the envelope mass is ~ 10~^Mq 
(Soker 1992). This implies that the post-AGB star loses 
~ few X 1O“^M0 during the relevant phase. From the 
discussion in §3.4 the mass loss rate is required to be 


~ fewx 1O“®M0 yr“^. According to the proposed model, 
therefore, the relevant phase lasts ~ 1000 yrs. From this 
discussion it seems quite plausible that the strongly in¬ 
teracting WD systems proposed in this paper comprise 
most, or even all, of the [O HI]-luminous PNs in galax¬ 
ies with an old stellar population. This might account 
for the finding that the cutoff luminosity of the PNLF in 
such galaxies is the same as that of galaxies with a much 
younger stellar population. 

4. SUMMARY 

Motivated by observations described in §2, I proposed 
that some of the [O HI]-luminous PNs are actually proto- 
PNs, where a WD companion accretes mass at a rela¬ 
tively high rate from the post-AGB star, sustains con¬ 
tinuous nuclear burning, and ionizes the dense nebula. 
Such systems are actually symbiotic type nebulae, sim¬ 
ilar to supersoft X-ray sources in many respects. How¬ 
ever, because of the dense nebula and large distances 
to the external galaxies, these are classified as PNs. In 
young stellar populations these strongly-interacting WD 
binaries account for a small fraction of [O HI]-luminous 
PNs, and single star models can account for the PNLF 
(e.g., Marigo et al. 2004). However, in old-stellar pop¬ 
ulations, such as in many elliptical galaxies, these bi¬ 
nary PNs might account for most, or even all, of the 
[O HI]-luminous PNs (§3.6). This could be the solution 
to the puzzling constant cutoff (maximum) [O HI] A5007 
luminosity of the PNLF across different galaxy types. 
The constant cutoff [O HI] A5007 luminosity is accounted 
for by the behavior of the accreting WDs: (i) Like for 
single central stars of PNs (Marigo et al. 2004) most 
WDs are concentrated around mass Mwd ~ O.6M0, 
with very small number having Mwd > 0.7Mq. (ii) 
For a given WD, the luminosity does not increase above 
a maximum value (van den Heuvel et al. 1992), which 
is Lwd = 1O^L0 for Mwd = 0-7Mq. (Hi) For higher 
luminosities the WD accretion rate should be higher, re¬ 
quiring a higher mass loss rate by the post-AGB star. 
This in turn implies that the ionized region is small and 
obscured. These three properties “conspires”, I propose, 
to give a constant cutoff [O HI] A5007 luminosity, some¬ 
what similar to the explanation of single stars models 
(Marigo et al. 2004). I note that the luminosity of these 
systems declines monotonically as the mass loss rate from 
the post-AGB star drops, hence accretion rate decreases. 
Later the post-AGB central star heats-up enough to ion¬ 
ize the nebula, and luminosity can rise again. These sys¬ 
tems, therefore, fill the entire bright part of the PNLF, 
and merge with the distribution of single stars. 

This model predicts that in old stellar populations the 
[O HI]-luminous PNs are bipolar, despite the low mass 
of their progenitors. However, because PNs are [O HI] 
luminous when they are young, hence small, it is very 
challenging to spatially resolve such PNs. As some frac¬ 
tion of [O HI]-luminous PNs in young stellar populations 
are also expected to be proto-PNs with accreting WDs, it 
will be interesting to search for such WDs inside [O HI]- 
luminous PNs in our galaxy. 

This research was supported in part by the Asher Fund 
for Space Research at the Technion. 
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